Local spin density in two-dimensional electron gas with hexagonal boundary 
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The intrinsic spin-Hall effect in hexagon- shaped samples is investigated. To take into account the spin-orbit 
couplings and to fit the hexagon edges, we derive the triangular version of the tight-binding model for the linear 
Rashba [Sov. Phys. Solid State 2, 1109 (I960)] and Dresselhaus [Phys. Rev. 100, 580 (1955)] [001] Hamil- 
tonians, which allow direct application of the Landauer-Keldysh non-equilibrium Green function formalism to 
calculating the local spin density within the hexagonal sample. Focusing on the out-of-plane component of 
spin, we obtain the geometry-dependent spin-Hall accumulation patterns, which are sensitive to not only the 
sample size, the spin-orbit coupling strength, the bias strength, but also the lead configurations. Contrary to the 
rectangular samples, the accumulation pattern can be very different in our hexagonal samples. Our present work 
provides a fundamental description of the geometry effect on the intrinsic spin-Hall effect, taking the hexagon as 
the specific case. Moreover, broken spin-Hall symmetry due to the coexistence of the Rashba and Dresselhaus 
couplings is also discussed. Upon exchanging the two coupling strengths, the accumulation pattern is reversed, 
confirming the earlier predicted sign change in spin-Hall conductivity. 

PACS numbers: 72.25.Dc, 73.23.-b, 71.70.Ej, 85.75.Nn 
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I. INTRODUCTION 



Spintronics, a science investigating the mechanism of elec- 
tron spins, has stimulated general theoretical interests in and 
attempts of industrial application on how to control electron 
spins. 1,2 By coupling the electron spin to the charge degrees of 
freedom, the spin-orbit (SO) interaction makes manipulating 
spin electrically possible. In particular, the spin-Hall effect 
(SHE), originating from such interaction, has attracted lots of 
researchers' attention since it induces a pure transverse spin 
current simply as a response to the longitudinal charge cur- 
rent. 

According to the sources of SO interaction, two types of 
the SHE can be identified: (i) extrinsic SHE^i^ with SO 
coupling being sourced by the external impurities, and (ii) 
intrinsic SHE (ISHE) with inherent SO coupling in the sys- 
tem. In both effects, the SO interaction subjects the spin-up 
(+z) and spin-down (—z) electrons respectively to opposite 
forces perpendicular to their transport direction. As a result, 
up-spin piles up on one side while down-spin on the other 
of the sample (in x-y plane), forming the so-called spin-Hall 
accumulation (SHA). Inspired by the recent breakthrough in 
experiments, 7 these two effects are intensively under investi- 
gation in bulk semiconductors, 8 ' 9 - 10 ' 11 - 12 - 13 ' 14 two-dimensional 
semiconductor heterostructureS j 15 i 16 ' 17 ' 18 ' 19 ' 20 ' 21 ' 22 - 23 ' 24 i 25 i 26 as 
well as the laterally confined quantum wiresj 2 - 7 ^ Moreover, 
due to the stronger magnitude of spin current in the intrin- 
sic case than in the extrinsic case (several orders greater), 
the ISHE has attracted intensive interests in theoretical study 
since the pioneering proposals*^ Nevertheless, most inves- 
tigations focus on the rectangular geometry. The affection of 
the boundary in other shapes is rarely discussed. Character- 



izing the SHE in finite system 
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this boundary de- 



gree of freedom as well as the varieties in bias configurations 
offer alternative ways for manipulating spin, and thus can fur- 
ther generate other applications. 



In this paper we study the geometry effect by considering 
a clean two-dimensional electron gas (2DEG) with hexagonal 
shape, which gives us more degrees of freedom in contacting 
the sample. Inside the 2DEG, two kinds of well-known and 
commonly referred SO couplings are considered: the Rashba 
SO (RSO) and the Dresselhaus SO (DSO) couplings. The for- 
mer is related to the inversion asymmetry of the structure 29 
with its coupling strength adjustable via the gate voltager^i 
the latter is caused by bulk inversion asymmetr y 32 ! 33 with 
coupling strength depending on the material. 34,35 Employing 
the Landauer-Keldysh (LK) Green function method in real 
space, 22 ' 23 ! 36 we analyze the z component SHA pattern inside 
the hexagon-shaped Landauer setup with two to six terminals, 
each of which may be contacted by a semi-infinite ideal lead. 
Electrochemical potential differences between these leads in- 
duce longitudinal charge current among terminals. Contrary 
to the SHE in rectangular samples with two head-to-tail leads, 
the SHA patterns in our hexagonal 2DEG preserve the funda- 
mental spin-Hall symmetry: the spins accumulate symmetri- 
cally (same magnitude) but oppositely (different sign) about 
the charge current flow direction, which is in general not 
straight in our hexagonal samples. Moreover, presented ac- 
cumulation patterns are found to depend on (i) SO interaction 
strength, (ii) bias strength, (iii) sample size, (iv) lead configu- 
ration, and (v) SO interaction type, among which we put more 
emphasis on the last three. For factors (iii) and (v), a SHA 
reversal effect (SHARE) is identified and discussed in detail. 
In particular, upon interchanging the RSO and DSO coupling 
strengths for factor (v) the whole SHA pattern is reversed, 
confirming the sign change in the spin-Hall conductivity as 
previously predicted in Ref. |20|. As for factor (iv), the ac- 
cumulation pattern in certain lead configuration may be very 
different from those in rectangular samples. 

This paper is organized as follows. In order to apply the 
real space tight-binding Hamiltonian in the LK formalism and 
fit the hexagon edges with appropriate discrete spatial points, 
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we derive the triangular lattice version for the linear Rashba 
and Dresselhaus tight-binding model in Sec. Q]] This version 
also enables further study of the geometry effect on the SHE 
with other shapes such as triangle, trapezoid, diamond, paral- 
lelogram, arrow, etc. Section [ill] gives a brief review on the 
LK formalism. Numerical results for SHA patterns based on 
the techniques introduced in Sec. [TT] and Sec. [HI] in various 
hexagonal samples will be reported and discussed in Sec. IIV1 
in the zero temperature limit. We finally conclude in Sec. [V] 



II. TIGHT-BINDING (FINITE DIFFERENCES) MODEL IN 
TRIANGULAR LATTICE WITH RSO AND DSO 
COUPLINGS 

In this section, we approach the original Hamiltonian by 
discretizing the continuous space into lattice-like points with 
the finite difference method. The following results are thus re- 
liable as the electron wave length is much longer than the dis- 
tance between nearest points. This method exactly treats the 
Hamiltonian when the realistic lattice points match the ones 
in finite differences, and it is also a very general approach for 
any lattice structures. With the constructed triangular lattice, 
we intend to answer the essential questions on how the sample 
boundary, which is basically controllable in today's technol- 
ogy, affects the spin accumulation. 

Consider a hexagon-shaped 2DEG (set in the x-y plane and 
grown along [001]) described by the Hamiltonian 



H = H K +H R + H D , 



(1) 



where H K — p 2 / '2m* (with m* being the electron effective 
mass) is the kinetic energy, and 



H R = — (p y <7 x - p x a y ) 



H 



(PxOx - PyOy 



(2) 

(3) 



denote the RSO and DSO Hamiltonians with coupling 
strengths a and (3, respectively. Here p x and p y are the 
two independent contributions to the two-dimensional mo- 
mentum operator: p — (p x ,p y ), and the Pauli matrices 
& = {vxiCyiCz) are used. To describe the system in real 
space, one can apply the tight-binding method to extract Ti 
in the real space lattice basis. However, to fit the hexagonal 
boundary, one has to discretize the space in triangular lattice 
structure, as shown in Fig. [T] where the lattice points are la- 
beled as 0, 1, • • • ,6, while A, B, ■ ■ ■ , G are auxiliary points. 
Using the finite difference method 36 and labeling the wave 
function at position r as ip T , the general expressions for first 
and second derivatives are approximated by 



— — 

Sr dr r 



'St d 
Sr dr 



Sr 

(Sr/2) 2 



(4) 



(5) 




FIG. 1: Schematic of the triangular lattice structure. On each lattice 
point, the electron has six nearest neighbors to hop. 



with r = (ri + Y2) /2 and <5r = ri — Y2- Thus the operation 
of these momentum operators can be expressed as 



Pxi>o = 

PylpO 



i h 

(Vh 

a 



-ih 



(1>i 



■ iPb) , 



(6a) 
(6b) 



with a being the lattice constant. Note the different denom- 
inators in Eqs. (|6ab and d6bl due to the different distances 
between HD and FB. In order to construct the matrix repre- 
sentation of the Hamiltonian in the real space triangular lattice 
basis, our goal is to express the operation of Tt on the wave 
function at point 0: ipo, in terms of those on its nearest neigh- 
bors: -01,^2, • • • Re- 
starting with the kinetic energy term, we take advantage 
of the isotropy of TL K , which implies that the six directions 
(from site to 1, 2, • • • , 6) in Ti K should take the same form. 
From Eq. ©, we therefore have two additional momentum 
operators p\ along EA and pi along GC with their operations 
on ijjQ given by 



pi ^0 
P21P0 



-ih 

a 

-ih 



(A 



-0c) • 



(7a) 
(7b) 



Combining the three equal contributions from p x ,Pi, and P2 
[by using Eqs. d6ab , ( |7al >, and jTbl l and defining the hopping 
energy to = h 2 /2m*a 2 , straightforward substitution and re- 
arrangement give 



P\ 



Pi 



2m* 

-t Q [{ih - 

("03 + 1p6 
-t [W>2 - 
("03 + ^6 



■00 

- v 4 - 2^0) + (ih +i>5- 2^0) 

- 2^0)] (8) 

- 0> 4 - 2lp D ) + ("01 + fa ~ 2Vi?) 

- 2?Ao)] - 2t (Vd +rl>B~ 2^o) • 



From the above equation, together with Eq. Q, we deduce 

[{p\+pl+p 2 x )/2m*}^ = (3/2)[(p»+pg)/2m*]Vo, 
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TABLE I: Matrix elements of the kinetic energy term TL K , the linear 



Rashba term TL R , and the linear Dresselhaus term Tt D in Eq. QJ in 
the triangular version of tight-binding model. 



and finally obtain 



hK = 2(pi+p 2 2 + P 2 x ) 



2m* 



(9) 



Note that here the long wave length limit p y (ipn + ipjj) /2» 
p 2 (tpo) is assumed. 

Next we seek for the matrix representation for 7i R and H D . 
Similar to the expression for H K , we need to express the op- 
erations p x ipo and Pyipo in terms of 0>i , 0>2, • • • , "06 ■ This can 
be done by noting 



p ^^b_+jPf_ _ _^/_ ^ _ ^ + ^ _ ^ _ 



(10) 



Py — = —j=- [("04 - 02) + (05 - 0l)] ~ PylpO, 



\/3a 



(11) 



where the long wave length limit is again assumed. Thus the 
operation on ipo due to terms required in the linear Rashba 
and Dresselhaus Hamiltonians given in Eqs. d2J and (0 can 
be completely expressed in terms of ipi , ip2, • • • ,ip6- 

Using Eq. © with Eq. © and Eqs. (fTOt— (ITTb. we ob- 
tain the finite difference representation of the Hamiltonian 

(/Hip) — E n =t'^(o,n)'/ , iT Since the matrix elements are 
non-vanishing only for nearest hopping: 

= ^(o,o)> if i = 3 

= W(o,n)) n G {1) 2, -.6}, if i, j are nearest neighbors 
= 0, otherwise, 



this expression is equivalent to the tight-binding model 
-^tbm _ j^. H(ij-)c\cj, with (c) being the creation (an- 
nihilation) operator. Defining the RSO and DSO hopping pa- 
rameters as t R = a /2a and t D = (3 /2a, respectively, the 
on-site energy Wmo) anc l hopping matrix elements H(o,n) ^ e 
summarized in Table U with each row corresponding to dif- 
ferent hopping direction referring to Fig. [TJ For example, 



the first row, labeled by (0, 0) represents the on-site energy 
or the energy 7~C(pm hopping from site to 0, while the sec- 
ond row, labeled by (0, 1), represents the energy W(o,i) hop- 
ping from site 1 to so that Table Q] gives Tiro o) = 4ip an d 



^(0,1) = 

it R <7 y /2) 



n/K 

(0,1) 



~ H (0, 



n/D 

" rt (0,l) 



-2t /3+(it R a x /V3- 



{it D a y /V3 - it D a x /2), etc. 



ffl. LANDAUER-KELDYSH FORMALISM 

In this section we give a brief review on the LK formal- 
ism, namely, the Keldysh non-equilibrium Green function 
formalism 37 ' 38 ' 39 ' 40 applied on Landauer multiterminal setups. 
The following review, which is in general applicable for meso- 
scopic systems with any shapes, is mainly based on Ref. |23|. 

Consider a hexagon-shaped Landauer setup, each side of 
which may be contacted by a semi-infinite two-dimensional 
ideal lead. In the following calculation, the leads are as- 
sumed to be in their own thermal equilibrium all the time 
(before and after contact), while the hexagonal conductor is 
in a non-equilibrium state, meaning that the electrons therein 
are not distributed according to a single Fermi-Dirac distri- 
bution. The occupation number of electron with spin a on 
site m at time t is determined by the diagonal matrix element 
(t, t; t s ) of the lesser Green function 



u mm,(T(T 



G 



m m.cr'a 



(t', t; t.) = I (cL, (i, t s ) CW (*', *»)) , (12) 



where cJ n(T (c mcr ) is the creation (annihilation) operator for the 
corresponding site m and spin index a. The lesser Green 
function given in Eq. (fT2l in general depends on the switch- 
ing time t s , at which the leads are brought into contact. (In 
the following expressions, t s will be set to — oo.) In the 
steady state, however, only the relative time variable t' — t = 
t is relevant, so that the lesser Green function can be ex- 
pressed in terms of its Fourier components G^, m , (r) = 

{2irhy 1 dEG^,^^ (E) e lET / h , leading to the occu- 
pation number 



{ c m c m) 



-G 



< 

m,m 



(r = 0) = — / dEG< m , aa (E) . 
2m J-oo 

(13) 

The solution to the lesser Green function matrix G < (E) is 
a kinetic equation G< (E) = G R {E) S< (E) G A (E) ,with 

G R (E) = [EI - HrffP 1 and G A (E) = [G R (E)] f being 
the retarded and advanced Green functions, respectively. The 
interactions between the conductor and the leads are included 
in the self-energy 



= [* 2 g P * (E)] 



{pi,Pi) 



(14) 



which yields the effective conductor Hamiltonian H e ff = H + 
J2 P ^p(E-eVp), where H is the conductor Hamiltonian with 
matrix elements Hu n summarized in TableU InEq. (fl4l >. the 
indices (i, j) on the left side and (pi.Pj) on the right side label 
the adjacent sites connecting the conductor and the lead p, and 
g^" (E) is the retarded Green function for the corresponding 
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isolated lead. The lesser self-energy S < , which generates the 
open channels, is given by 



N = 20 



E< (E) = -J2 [V P (E - eV p ) - S], (E - eV p )] 
p 

xf(E- eV p ) , 



(15) 



where / (E — eV p ) is the Fermi-Dirac distribution and V p is 
the applied bias voltage on lead p. Note that Eq. ([T5T l im- 
plies that G< is an anti-Hermitian matrix ensuring that the 
occupation number is a real quantity, while the non-Hermitian 
Hamiltonian H e ff yields the finite eigenstate lifetime of quasi- 
particle. 



IV. NUMERICAL RESULTS 

In this section we present numerical results for the z com- 
ponent of the local spin density (S z ) inside the hexagonal 
sample, made of InGaAs/InAlAs heterostuctures grown along 
[001]. Typical parameters: the electron effective mass m* — 
0.05m e (m e is the electron rest mass) and the lattice constant 
a = 3 nm (yielding the hopping energy to — 84.68 meV), 
will be taken, 30 and the Fermi energy Ep = — 3.8io close to 
the band bottom Ei, = — 4fo will be chosen (so that the long 
wave length approximation is valid). The size of the regu- 
lar hexagonal samples will range from N = 8 to N = 20, 
N being the number of sites per hexagon edge, while the 
applied bias on the leads will be set either eV p = +eV/2, 
eV p = —eV/2, or eV p — 0, to be denoted shortly as "+", 
"— ", and "0" in each presented figure, respectively. Note that 
here the electron is negatively charged as e = — |e|, so that 
eV p > implies V p < 0, meaning that electrons will flow 
from "+" to "— ", in all the presented figures. We will take 
eV = 0.4io and eV = 10 -3 <o, which will be referred to as 
high bias and low bias, respectively. Discussion with Rashba- 
type 2DEGs, for which t R = O. ltp an d t D = are set, will 
be first addressed. Not until Sec. II V CI will we change t R and 
turn on t D . Note that in most of our numerical results, we put 
parameters identical with Ref. 22, so that direct comparison 
and correspondence can be clearly identified. 



A. SHA: General feature 




x 10" 



x 10 ' 



FIG. 2: (Color online) Non-equilibrium spin-Hall accumulation (S z ) 
in a (a) = 20 sample under high bias, and a (b) N — 8 sample 
under low bias, (c) is the same as (b), except the reversed bias. The 
color bars calibrate (S z ) in units of h/2, and ± signs indicate the 
applied bias ±eV/2 on the corresponding lead. 



force 24 depends on not only the sign of S z , but also the system 
size relative to the Rashba spin precession length^ 



irh 2 irato 
2m* a ~ 2t R ' 



(16) 



which is the distance required by the RSO interaction to rotate 
the spin by it. With t R = O.lto here, Eq. ( TToT l yields the 
precession length Lso — 5ira w 15.7a. We will address this 
accumulation reversal effect a bit further later. Upon reversing 
the bias, Fig. |2|c) shows a flipped pattern, as expected. In 
short, Fig. [2] can be regarded as the hexagonal version of Fig. 
1 shown in Ref. 22, and the general feature of the intrinsic 
SHE is preserved. 



As mentioned in Sec. U the SHA pattern in hexagon-shaped 
samples can be very different from those in rectangular ones. 
Beginning with Fig. Ufa), we plot (S z ) in a TV = 20 sample 
with two head-to-tail leads under high bias. The accumulation 
pattern therein clearly shows the geometry effect: the edge 
peaks of the SHA follow the shape of the sample. Moreover, 
due to the head-to-tail lead configuration, the accumulation 
pattern is perfectly spin-Hall-symmetric about the central axis 
connecting the two leads: spin accumulation on one side is the 
mirror (with different sign) of that on the other. 

When reducing the sample size to N = 8, one can see 
that the accumulation pattern changes sign, as shown in Fig. 
|2jb), where low bias is applied. This is because the spin-orbit 



B. SHA: Symmetry investigation 

From now on we focus on the low-bias regime and, in this 
subsection, we put emphasis on different lead configurations. 
Various N = 12 samples with 2 leads, 4 leads, and 6 leads 
with different configurations will be shown. 

Beginning with the 6-lead cases, Figs. [3ja)-[3]c) show SHA 
patterns under a series of different biasing. For the non-head- 
to-tail biasing in Figs. [3ja) and Ob), one can see a curved 
spin-Hall symmetry which roughly follows the charge current 
flow. Moreover, such a bent SHA may induce (S z ) on edges 
with not necessarily the opposite signs, e.g., see the dark spots 
in the left-top and right-bottom corners shown in Fig. fSa). 
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FIG. 3: (Color online) Spin-Hall accumulation in N = 12 hexagonal 
samples with (a)-(c) 6 leads, (d)-(f) 2 leads, and (g)-(h) 4 leads, 
under low bias. In addition to the ± signs indicating the applied 
±eV/2 bias, the zero signs in (a)-(c) label the unbiased leads. Color 
bars calibrate (S z ) in units of (h/2) x 10~ 5 . 



For the head-to-tail-biased sample of Fig. [3| C X the straight 
spin-Hall symmetry is recovered, as expected. As a direct cor- 
respondence, we unplug those unbiased electrodes and show 
the accumulation patterns in Figs. [Ud)-(f). Similar proper- 
ties mentioned above are still valid, but each pattern becomes 
totally different, implying that the spin accumulation is sensi- 
tively affected by the configuration of attached leads. 

We now focus on Figs. [3jc) and 0T) and discuss the 
SHARE for these two distinct cases: open vs. closed bound- 
aries. For the open case Fig. 0c), one can clearly observe that 
the accumulation pattern reverses with different transverse 
length, which ranges from (N — 1) a = 11a (the hexagon 
edge length) to 2 (N — 1) a = 22a (the maximal length of the 
hexagon) here and covers the precession length Lso = 15.7a. 
Near the biased top and bottom electrodes, the spin accumula- 
tion signs can be well explained by the semiclassical SO force 
(to be discussed further later) proportional to (p x e 2 ) ® a z 
[see Eq. ( fT8l ) in Sec. IIVCI , which predicts right deflection 



for cr z = +1 and left deflection for a z = —1. At central 
regions where the transverse length exceeds Lso, the Rashba 
spin precession turns a z = ±1 into a z = =pl and thus flips 
the pattern. 

For the closed case Fig. [3jf), the SHARE is slightly dif- 
ferent, as one can see that the bright and dark regions do not 
swap locally. A basic difference between the open and closed 
cases is that the electron does not have the outflow degree of 
freedom in the latter case. Therefore, the criterion of whether 
the pattern is flipped or not lies on the maximum width of 
the whole sample, which amounts to 22a for this N = 12 
case and exceeds Lso already. Accordingly, the SO force pre- 
dicts right (left) deflection for a z = — 1 (a z = +1) elec- 
trons, in agreement with Fig. Of). Indeed, the critical size for 
the closed sample (with two head-to-tail leads) to exhibit the 
flipped accumulation patterns can be shown as N — 9, but we 
do not further show here. 

Figures [3jg) and [3jh) show accumulation patterns in the 
hexagon sample attached by 4 leads with two different bias 
configurations. As in the former case longitudinal bias is ar- 
ranged similar to the rectangular sample attached to two leads, 
regular spin-Hall accumulation is observed in Fig. [3jg)- In- 
terestingly, when we rearrange the applied bias as that in Fig. 
Oh), "local" SHA is obtained: spin accumulation in the whole 
sample is divided into two parts, subject to the two separate 
pairs of positive-negative electrodes. The reason to obtain 
such a divided accumulation pattern can be intuitively under- 
stood as that each pair of the electrodes are closely connected, 
so that electron spins flow directly from "+" to "— " at one 
side and are less affected by the electrode pair at the other 
side. 

A last remarkable point for Fig. [3] is the order of magni- 
tude of (S z ). Strictly speaking, one should refer the local spin 
density (S z ) to the spin accumulation only for those patterns 
in closed samples. When attaching unbiased transverse leads, 
electrons are free to leak out, such that the local spin density 
shows merely the spatial spin distribution. It is only when the 
leads at the transverse sides are unplugged that those spins 
deflected by the SO force accumulate at the transverse edges 
and that the local spin density showing obvious peaks near the 
edges deserve the name spin-Hall accumulation. In addition, 
such a leakage of the spins due to plugged leads (not neces- 
sary unbiased) may lower the order of magnitude of (S z ) as 
can be clearly seen by comparing the three groups in Figs. [3] 
(a)-(c), (d)-(f), and (g)-(h). 

C. Competition forces in-between Rashba and Dresselhaus 
interactions 

The semiclassical force provides a simple description of the 
SHE. Using the Heisenberg equation of motion with Hamil- 
tonian H given by Eqs. ([l])-©, the velocity can be obtained 

as 

v = r=i[r,fl] 

in 

= -^7 + r [{-OiUy + (3a x ) e x + (aa x - j3a y ) e y ] , (17) 
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leading to the SO-coupling-induced force 

^ * . m * r n 2m* , , , , 

F so = m*v =— [v, H] = -p- (a 2 - (p x ej a z . 

(18) 

Apart from the constant prefactor 2m* /ft', three transpar- 
ent properties of Fso given by Eq. ( fT8l can be identified. 
First, spin-up electrons (a z = +1) and spin-down electrons 
(<j 2 = — 1) feel opposite forces. Second, since Fso is al- 
ways perpendicular to the electron transport direction [Fso oc 
(p x e z ) _L p], no work is done by this force. In other words, 
the spin current caused by Fso is dissipationless. Third, the 
Rashba and Dresselhaus interactions contribute oppositely to 
Fso, such that they compete against each other. Thus the sign 
of Fso is also governed by (a 2 — /3 2 ), implying a reverse of 
the entire accumulation pattern when interchanging a with (3. 

To demonstrate this SHARE due to the exchange of 
(a,/3) ^ (/3,a), or (t R ,t D ) <-> (t D ,t R ), let us con- 
sider again the low-biased N = 12 hexagonal samples with 
two head-to-tail leads. We show SHA patterns for various 
(t R , t D ^j in Fig. |4] where the left and right columns are 
related by interchanging t R and t , and the SHARE in- 
duced by the DSO coupling is clearly seen: patterns with 
(t R , t j <-> (t D , £ fl ) are reversal of each other. One can also 
observe that whereas we put t R > t D in the left column and 
t R < t D in the right column, the reversed patterns mean the 
exchange between the spin-up and spin-down states, imply- 
ing a sign change in the spin-Hall conductivity as previously 
predicted^ 

Another interesting point is that the coexistence of the RSO 
and DSO couplings breaks the spin-Hall symmetry. When 
one of t R or t D dominates, such as Figs. HJa) orHfd), the 
SHA on one side is (nearly) an mirror (with different sign) of 
that on the other side; while the magnitude of t R and t D be- 
comes competitive, such as Figs. |4jc) or|4jf), the spin-Hall 
symmetry is broken, or distorted. Interestingly, this distortion 
effect does not simply kill the SHA pattern. Instead, the spin 
accumulation on only one side becomes vague while that on 
the other side becomes even stronger, as shown in Figs. Hfc) 
orSJf). However, for systems with t R w t D vanishing (S z ) 
inside the 2DEG is obtained (not shown here), as already im- 
plied in Eq. (H). 



V. CONCLUDING REMARKS 

In conclusion, the intrinsic SHE in hexagon-shaped 2DEG 
is investigated. The tight-binding model in triangular lattice 
which suits the hexagon edges is constructed (within the long 
wave length limit) for the linear Rashba and Dresselhaus [00 1 ] 
models. Applying the derived Hamiltonian (with matrix ele- 
ments summarized in Table HJ to the LK formalism, we have 
shown the SHA patterns in hexagonal samples under various 
conditions, including SO interaction strength, bias strength, 
sample size, lead configuration, and SO interaction types. The 
spin-Hall symmetry is identified and discussed in detail, while 
the reversal effect of the SHA pattern, which may be induced 
by the sample size and the competition between the RSO and 
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FIG. 4: (Color online) Spin-Hall accumulation (S z ) in units of 
(ft/2) x 10 for N = 12 hexagonal samples, each with two head- 
to-tail leads under low bias, for various t R 's and t D 's. Left column 
(a)-(c) and right column (d)-(f) are related by interchanging t R and 



DSO couplings, are discussed. Contrary to the SHA in rect- 
angular samples, we demonstrate here that the accumulation 
pattern can be very different due to the geometry effect. 

Before closing, three remarks are worthy of mention here. 
First we highlight again the main difference between the 
hexagonal and the rectangular samples. Although bending the 
spin-Hall symmetry by applying non-straight biasing can be 
done in the rectangular case, transport angles with multiples 
of 90 degrees lowered to 60 enlarge the option of lead con- 
figuration. Explicitly, the total number of ways of biasing the 
sample, e.g., for the two-lead case, is increased from C| = 6 
to Cf = 15, making more detailed analysis possible. (Note 
that in general, samples should be treated as crystallography- 
dependent, so that each biasing direction is nontrivial.) In ad- 
dition, the oblique side boundary varying with longitudinal 
position making the SHARE observable [Fig. |3]c)]. Such a 
reversal due to spin precession cannot be observed in one sin- 
gle rectangular sample. 

Secondly, we remark on the sample size dependence. There 
are two characteristic lengths relavant in our hexagonal sam- 
ples: spin precession length Lso and Fermi wave length Xp. 
As we have discussed the spin-orbit force, the former is the 
sign criterion of the spin accumulation pattern. As long as 
the device size is larger than Lso, the direction of the spin- 
orbit force, and hence the sign of (S z ) accumulating at the 
two sides, oscillate with the increase of the sample size. The 
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latter, Xp , modulates the spatial distribution of (S z ) and thus 
induces additional peaks similar to the standing waves. This 
means that in samples with the dimension of the order of or 
larger than Xp, such waves can be observed. In our calcu- 
lations, we have Ep — Ef, = 0.2io ~ 17 meV and thus 
Xf = 2irh/^/2m*Ep sa 13a, which is slightly shorter than 
Lso- A good example is to compare Fig. |2jb) and Fig. [3jf). 
Under exactly the same conditions (lead configuration, bias 
strength, etc.), the former with N — 8 shows no additional 
peaks, while in the latter with N = 12 the (S z ) distribution 
exhibits two peaks at each side due to wave function modula- 
tion. 

Finally, we stress again the nature of our calculation. It 
is the hexagonal boundary that makes discretizing the two- 
dimensional space inside the sample into a triangular lattice 
necessary. In the long wave length limit, where the Fermi 
energy is close to the band bottom so that crystal structure 



information is less important, our calculation simply reveals 
the free electron gas behavior. It is only when the Fermi level 
no longer approaches the band bottom that the band structure 
effect enters and the realistic crystal structure needs be taken. 

Our present work provides a fundamental description of the 
geometry effect on the intrinsic SHE, taking the hexagon as 
the specific case. Additional degrees of freedom on configur- 
ing the attached leads may imply alternative ways of manipu- 
lating spins and even other spin-Hall experiment setups. 
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